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Abstract 

We report the results from a study of a partial sample of ~ 2.3 X lO'^ 
TT^TT^TT^ decays recorded by the NA48/2 experiment at the CERN SPS, showing an 
anomaly in the tt^tt^ invariant mass (Mqo) distribution in the region around Moo = 
2m_|_, where m_|_ is the charged pion mass. This anomaly, never observed in previous 
experiments, can be interpreted as an effect due mainly to the final state charge 
exchange scattering process vr^vr^ vr'^vr" in — > 7r^7r"*"7r decay pp. It provides 
a precise determination of aQ—a2, the difference between the vrvr scattering lengths in 
the isospin 7 = and 1 = 2 states. A best fit to a rescattering model corrected for 
isospin symmetry breaking gives {ao — a2)m^ = 0.268±0.010 (stat.) ±0.004 (syst.), 
with additional external uncertainties of ±0.013 from branching ratio and theoretical 
uncertainties. If the correlation between oq and 02 predicted by chiral symmetry 
is taken into account, this result becomes (oq — a2)m^ = 0.264 ± 0.006 {stat.) ± 
0.004 (syst.) ±0.013 (ext.). 

To be published in Physics Letters B 
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1 Introduction 



The NA48/2 experiment at tlie CERN SPS is searching for direct CP violation in K"^ 
decay to three pions. The experiment uses simultaneous and K~ beams with a 
momentum of 60 GeV/c propagating along the same beam line. Data have been collected 
in 2003-04, providing samples of ~ 4 x 10^ fully reconstructed K"^ -k^-k^-k^ and ~ 10^ 
TT^vr'^Tr^ decays. Here we report the results from a study of a partial sample 
of ~ 2.3 X 10'' K"^ — > TT^vr'^Tr^ decays recorded in 2003, showing an anomaly in the 
7r°7r° invariant mass (Mqo) distribution in the region around Mqo = 2m+, where m+ is 
the charged pion mass. This anomaly, never observed in previous experiments, can be 
interpreted as an effect due mainly to the final state charge exchange scattering process 
TT+TT" — > Ti^TT^ in 7[^7i~^TT~ dccay pp. A best fit to a rescattering model provides 

a precise determination of aQ — a2, the difference between the S-wave vrvr scattering lengths 
in the isospin 7 = and 1 = 2 states. 

2 Beam and detectors 

The two simultaneous beams are produced by 400 GeV protons impinging on a 40 cm 
long Be target. Particles of opposite charge with a central momentum of 60 GeV/c and 
a momentum band of ±3.8% produced at zero angle are selected by a system of dipole 
magnets forming an "achromat" with null total deflection, focusing quadrupoles, muon 
sweepers and collimators. With 7 x 10^^ protons per burst of ~ 4.5 s duration incident 
on the target the positive (negative) beam flux at the entrance of the decay volume is 
3.8 X 10^ (2.6 X 10^) particles per pulse, of which ~ 5.7% (~ 4.9%) are K+ {K'). The 
decay volume is a 114 m long vacuum tank with a diameter of 1.92 m for the first 66 m, 
and 2.4 m for the rest. 

Charged particles from K"^ decays are measured by a magnetic spectrometer consisting 
of four drift chambers [S] and a large-aperture dipole magnet located between the second 
and third chamber. Each chamber has eight planes of sense wires, two horizontal, two 
vertical and two along each of two orthogonal 45° directions. The spectrometer is located 
in a tank filled with helium at atmospheric pressure and separated from the decay volume 
by a thin (0.0031 radiation lengths, Xq) Kevlar window. A 16 cm diameter vacuum tube 
centered on the beam axis runs the length of the spectrometer through central holes in 
the Kevlar window, drift chambers and calorimeters. Charged particles are magnetically 
deflected in the horizontal plane by an angle corresponding to a transverse momentum kick 
of 120MeV/c. The momentum resolution of the spectrometer is (t(p)/p = 1.02%©0.044%]? 
(p in GeV/c), as derived form the known properties of the spectrometer and checked 
with the measured invariant mass resolution of tt^tt^tt^ decays. The magnetic 

spectrometer is followed by a scintillator hodoscope consisting of two planes segmented 
into horizontal and vertical strips and arranged in four quadrants. 

A liquid Krypton calorimeter (LKr) j4j is used to reconstruct 7r° 77 decays. It is 
an almost homogeneous ionization chamber with an active volume of ~ 10 iv? of liquid 
krypton, segmented transversally into 13248 2 cm x 2 cm projective cells by a system 
of Cu-Be ribbon electrodes, and with no longitudinal segmentation. The calorimeter is 
27 Xq thick and has an energy resolution a{E)/E = 0.032//E © 0.09/ E © 0.0042 (E 
in GeV). The space resolution for single electromagnetic shower can be parametrized as 
c"x = o'y = 0A2/\fE © 0.06 cm for each transverse coordinate x, y. 



A neutral hodoscope consisting of a plane of scintillating fibers is installed in the LKr 
calorimeter at a depth of ~ 9.5 Xq. It is divided into four quadrants, each consisting of 
eight bundles of vertical fibers optically connected to photomultiplier tubes. 

3 Event selection and reconstruction 

The —* vr^vr^vr^ decays are selected by a two level trigger. The first level requires 
a signal in at least one quadrant of the scintillator hodoscope in coincidence with the 
presence of energy depositions in LKr consistent with at least two photons. At the second 
level, a fast on-line processor receiving the drift chamber information reconstructs the 
momentum of charged particles and calculates the missing mass under the assumption 
that the particle is a tt^ originating from the decay of a 60 GeV/c travelling along 
the nominal beam axis. The requirement that the missing mass is not consistent with 
the vr" mass rejects most of the main — > tt^tt^ background. The typical rate of this 
trigger is ~ 15, 000 per burst. 

Events with at least one charged particle track having a momentum above 5 GeV/c, 
and at least four energy clusters in LKr, each consistent with a photon and above an energy 
threshold of 3 GeV, are selected for further analysis. In addition, the relative track and 
photon timings must be consistent with the same event within the experimental resolution 
(~ 1.5ns). The distance between any two photons in LKr is required to be larger than 
10 cm, and the distance between each photon and the impact point of any track on LKr 
must exceed 15 cm. Fiducial cuts on the distance of each photon from the LKr edges and 
centre are also applied in order to ensure full containment of the electromagnetic showers 
and to remove effects from the beam pipe. Finally, the distance between the charged 
particle track and the beam axis at the first drift chamber is required to be larger than 
12 cm. 

At the following step of the analysis we check the consistency of the surviving events 
with the — ^ vr^vr^vr^ decay hypothesis. We assume that each possible pair of photons 
originates from 7r° — 77 decay and we calculate the distance Dik between the 7r° decay 
vertex and the LKr: 

_ ^jEiEk[{xi - XkY + {Vi - VkY] 

Dik — 

mo 

where Ei,Ek are the energies of the i-th and k-th photon, respectively, Xi, yi, Xk, Vk are the 
coordinates of the impact point on LKr, and mo is the 7r° mass. Among all photon pairs, 
the two with the smallest Dik difference are selected as the best combination consistent 
with the two mesons from tt^tt^tt^ decay, and the distance of the decay 

vertex from the LKr is taken as the arithmetic average of the two Dik values (it can be 
demonstrated that this choice gives the best ir^n^ invariant mass resolution near thresh- 
old). Fig. ^ shows the invariant mass distribution of the system consisting of the two 7r° 
and a reconstructed charged particle track, assumed to be a n"^. This distribution is dom- 
inated by the peak, as expected. The non Gaussian tails originate from unidentified 
TT^ fi^ in flight or wrong photon pairing. The final event selection requires that the 
7r^7r°7r° invariant mass differs form the mass by at most ±6 MeV. This requirement 
is satisfied by 2.287 x 10^ events. The fraction of events with wrong photon pairing in 
this sample is ~ 0.25%, as estimated by a high-statistics fast Monte Carlo simulation of 
—>■ 7r^7r°7r° decays which takes into account the momentum distribution of the three 
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Figure 1: Invariant mass distribution of reconstructed tt^tt^tt^ candidate events. The arrows indicate 

ttie selected mass interval. 

pions, 7r° 77 decay kinematics and the effect of the detector acceptance and resolution. 



4 Cusp anomaly in the tt^tt" invariant mass distribution 

Figure El shows the distribution of the square of the tt^tt^ invariant mass, Mqq, for the 
final event sample. This distribution is displayed with a bin width of 0.00015 (GeV/c^)^, 
with the 51** bin centered at Mqq = (2m+)^ (as discussed below, the bin width is chosen 
to be smaller than the Mqq resolution). A sudden change of slope near Mqq = (2m+)^ = 
0.07792 (GeV/c^)^ is clearly visible. Such an anomaly has not been observed in previous 
experiments. 

The Dalitz plot distribution for —>■ Tr'^Tr^Tr" decays is usually parametrized by a 
series expansion in the Lorentz-invariant variable u = [s^ — SQ)/m^, where Sj = {Pk — 
PiY (i=l,2,3). So = (si + S2 + Ss)/3, Pk (Pi) is the K{tt) four-momentum, and i = 
3 corresponds to the vr^ [5^. In our case S3 = Mqq, and Sq = (m|^ + 2mQ + m^)/3. 
We have used this parametrization in a fast Monte Carlo simulation of K"^ — >■ vr^Tr^Tr" 
decays with the same detector parameters used in previous NA48 analyses jS]. This 
simulation takes into account most detector effects, including the trigger efficiency and 
the presence of a small number (< 1%) of "dead" LKr cells. For any given value of the 
generated vr^vr^ invariant mass the simulation provides the detection probability and the 
distribution function for the reconstructed value of Mqq. This allows the transformation of 
any theoretical distribution into an expected distribution which can be compared directly 
with the measured one. 




X 10^ 




0.08 0.09 0.1 



0.11 0.12 0.13 

(GeV /cy 



Figure 2: Distribution of Mqq, the square of the tt^tt'^ invariant mass. The insert is an enlargement of a 
narrow region centered at Mqq = (2to+)^ (this point is indicated by the arrow). The statistical error bars 

are also shown in these plots. 

Figure shows the expected Mqq resolution (r.m.s.) as a function of Mqq, together 
with examples of Mqq distributions for five generated Mqq values. The Mqo resolution 
is the best at small Mqq values, varying between ~ 0.4 MeV/c^ near Mqq = 2mQ, and 
~ 1.4 MeV/c^ at the end of the Mqq allowed range. It is 0.56 MeV/c^ at Mqq = 2m+. 
A plot of the overall detector acceptance as a function of the generated Mqq value, as 
predicted by the Monte Carlo simulation (see Fig. IHb), shows no structure in the Mqq 
region where the sudden change of slope is observed in the data. 

We have tried to fit the distribution of Fig. |21 in the interval 0.074 < Mqq < 
0.097 (GeV/c^)^ using the distribution predicted by the Monte Carlo simulation with 
a matrix element as given in ref. P]: 

A^o = 1 + ^gou (1) 

In this fit the free parameters are go and an overall normalization constant. Because of 
the anomaly at Mqq = (2m+)^, it is impossible to find a reasonable fit to the distribution 
of Fig. 121 (the best fit gives = 9225 for 149 degrees of freedom). However, fits with 
acceptable values are obtained if the lower edge of the fit interval is raised few bins above 
Mq2q = (2m+)2. As an example, a fit in the interval 0.07994 < Mq^q < 0.097 (GeV/c^)^, 
with the lower edge only 0.002 (GeV/ c^)^ above (2m+)^, gives = 133.6 for 110 degrees of 
freedom. This fit gives Qq = 0.683 ±0.001 (statistical error only), in reasonable agreement 
with the present world average, go = 0.638 ± 0.020 [5j (it should be noted, however, 
that the matrix element used here has not the same form as that used in ref. j3]). 
The quahty of this fit is illustrated in Fig. HI which displays the quantity A = (data - 
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Figure 3: a) Expected Mqq resolution (r.m.s. in MeVje?) versus generated Mqq (full line histogram), 
together with Mqq distributions for five generated values of Mqq; b) Acceptance versus Mqq (see text). 

The point Mqq = (2m+)^ is indicated by the arrow. 
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Figure 4: A = (data - fit) /data versus Mqq. The point Mqq = (2rn+)^ is indicated by the arrow. Also 

shown is the Mgg region used in the fit. 

fit)/data as a function of M^g for the fit region 0.07994 < M^g < 0.097 {GeV/c^ f and 
also for Mqq < 0.07994 (GeV/c^)^, where the prediction with the same parameters is 
extrapolated. 

Figure m shows that, in the region Mqq < (2m+)^, the data fall below the prediction 
based on the same parameters obtained from the fit region. The total number of events 
in the first 50 bins of the data is 7.261 x 10^, while the extrapolated prediction gives 
8.359 X 10^ events. 

In order to investigate the origin of this "deficit" of events in the data we have studied 
the event shape distributions in two 20 bins wide intervals, one just below and the other 
just above Mqq = (2m+)^. Since Mqq is computed using only information from the LKr 
calorimeter, we consider only photon cluster parameters. We denote the distributions of 
measured photon energy and distances in these two intervals as /_ and J+, respectively, 
and we compare the J+ / /_ ratios with those predicted by the simulation after normalizing 
/_ and /+ to the same area. 

These ratios (see Fig. and Ej) show that the shapes of all distributions for the two 
Mqq intervals, as measured in the data, are in excellent agreement with the Monte Carlo 
predictions. In addition, no difference is observed between and K~ nor between the 
data taken with opposite direction of the spectrometer magnetic field. The simulation 
also shows that the Mqq distribution of the events affected by wrong photon pairing has no 
local structures over the whole Mqq range. We conclude that the Monte Carlo simulation 
describes correctly the Mqq dependence of the detection efficiency in the region around 
Mqq = (2m+)^, and the "deficit" of events in the data in the region Mqq < (2m+)^ is due 
to a real physical effect. 
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Figure 5: Data (points with error bars) - Monte Carlo (histogram) comparison of the ratio of normalized 
photon energy distributions between events with Mqq > (2m+)^ and Mqq < (2m+)^ (see text) 




Figure 6: Data (points with error bars) - Monte Carlo (histogram) comparison of ratios I+/I- of 
normalized distance distributions between events with Mqq > (2m+)^ and Mqq < (2m+)^ (see text), a) 
min. r^: distance (cm) between LKr centre and closest photon; b) max. r^: distance (cm) between LKr 
centre and farthest photon; c) min. d^^: minimum distance (cm) between photons at LKr; d) d-y-track- 
minimum distance (cm) between photons and tracks at LKr. 



5 Interpretation and determination of the tttt scattering lengths 



The sudden change of slope observed in the Mqq distribution at Mqq = (2m+)^ (see Fig. 
E} suggests the presence of a threshold "cusp" effect from the decay K"^ — > Tr^'^Tr^vr" 
contributing to the K"^ -k^-k^-k^ amplitude through the charge exchange reaction 
TT+TT^ tt'^tt". The presence of a cusp at Mqq = (2m+)^ in tt^tt" elastic scattering 
due to the effect of virtual tt+tt" loops has been discussed first by Meissner et al. [7j. For 
the case of —>■ tt^tt^tt'^ decay Cabibbo has proposed a simple rescattering model pQ 
describing the K"^ —>■ 7r='=7r°7r° decay amplitude as the sum of two terms: 



M{K^ 7r±7r%") =Mo + Mi 

where Aio is the "unperturbed amplitude" of Eq. ^ and A^i is the contribution from 
the TT^n^n^ decay amplitude through 7i~^7r~ 7r°7r° charge exchange, with the 

renormalization condition Aii = at Mqq = (2m_(_)^. The contribution Aii is given by 



Ml = -2a^m+M- 



M, 



00 



(2) 



where is the S-wave tc^tt^ charge exchange scattering length (threshold amplitude), 
and is the known tt^tt+tt^ decay amplitude at Mqq = 2m+. Aii changes from 

real to imaginary at Mqq = 2m+ with the consequence that Aii interferes destructively 
with M.Q in the region Mqq < 2m+, while it adds quadratically above it. In the limit of 
exact isospin symmetry Qx = {qq — a2)/S, where qq and 02 are the S-wave vrvr scattering 
lengths in the 1 = and 1 = 2 states, respectively. 

In this simple rescattering model there is only one additional parameter, axfn^. A 
fit to the Mq2q distribution in the interval 0.074 < Mq^q < 0.097 (GeV/c^)^ using axm+ 
as a free parameter gives = 420.1 for 148 degrees of freedom. The quality of this fit 
is illustrated in Fig. [7^ which displays the quantity A defined in section 4 as a function 
of Mqq. One can see that this model provides a much better but still unsatisfactory 
description of the data. In particular, the data points are systematically above the fit in 
the region near Mqq = (2m+)^. 

Recently, Cabibbo and Isidori j2] have proposed a more complete formulation of the 
model which takes into account all rescattering processes at the one-loop and two-loop 
level. In this formulation the matrix element for — >■ ir^n^ir^ decay includes several 
additional terms which depend on five S-wave scattering lengths, denoted by a^, a++, 
a+_, a+O; o-oo, and describing tt+tt^ tt'^tt'^, 7r+7r+ 7r+7r+, tt+tt^ tt+tt^, tt+tt^ — >■ 
TT+TT^, 7r°7r° Ti^TT^ Scattering, respectively. In the limit of exact isospin symmetry these 
scattering lengths can all be expressed as linear combinations of Qq and a2- 

At tree level, omitting one-photon exchange diagrams, isospin symmetry breaking 
contributions to the elastic tttt scattering amplitude can be expressed as a function of 
one parameter e = (m^ — mQ)/m^ = 0.065 |H|. In particular, the ratio between the 
threshold amplitudes a^, cl++, a+o, aoo and the corresponding isospin symmetric ones 
- evaluated at the vr^ mass - is equal to 1 — e for tt'^tt'^ 7T~^7r~^, n'^'n'^ — » tt'^tt'^, vr^vr^ — ^ 
tt'^tt^; 1 + e for Tr"*"?!" — n~^7r~] and 1 + e/3 for 7r"'"7r~ —>■ tt^tt^. These corrections have 
been applied to the rescattering model of ref. j2] in order to extract oq and 02 from the 
fit to the data. 



In the model of ref. [2] the matrix element for 7r^7r°7r° decay includes terms 

which depend on both independent kinematic variables (Mqo and M_|_o, the invariant mass 
of the 7r^7r° pair) requiring, therefore, a fit to the two-dimensional Dalitz plot. We have 
performed an approximate fit to this model by calculating these terms at the average value 
of M^Q for each value of Mqq. This fit has five free parameters: (ao — a2)m+, a2m+,go, 
a quadratic term of the form 0.5h'u^ added in equation ^ and an overall normalization 
constant. The quality of the fit (x^ = 154.8 for 146 degrees of freedom) is shown in 
Fig. [Tjp. A better fit (x^ = 149.1 for 145 degrees of freedom, see Fig. [7\:) is obtained by 
adding to the model a term describing the expected formation of tt+tt" atoms ( "pionium" ) 
decaying to vr'^vr^ at Mqo = 2m+. The best fit value for the rate of vr^+pionium 
decay, normalized to the tt^tt+tt^ decay rate, is (1.61 ±0.66) x 10~^, in reasonable 

agreement with the predicted value ~ 0.8 x 10~^ jH]. 

The rescattering model of ref. j2] does not include radiative corrections, which are 
particularly important near Mqo = 2m+, and contribute to the formation of tt'^7t~ atoms. 
For this reason we prefer to exclude from the final fit a group of seven consecutive bins 
centered at Mqo = 2m_|_. The quality of this fit (x^ = 145.5 for 139 degrees of freedom) 
is illustrated in Fig. [7|i, which shows the small excess of events from pionium formation 
in the bins excluded from the fit. Table ^ lists the best fit values of the parameters, 
as obtained by two independent analyses which use different event selection criteria and 
different Monte Carlo simulations to take into account acceptance and resolution effects 
(the analysis described so far is denoted as Analysis A; analysis B uses a simulation of the 
detector based on GEANT jTHI). We take the arithmetic average of these values as the 
measurement of these parameters, and one half of the difference between the two values 
as a systematic uncertainty from acceptance calculations. In both analyses changing the 
selection criteria never leads to variations of the best fit parameters larger than these 
uncertainties. 



Parameter 


Analysis A 


Analysis B 


Arithmetic average 


(oo - a2)m+ 


0.269 ±0.010 


0.268 ±0.010 


0.268 ±0.010 


a2m^ 


-0.053 ±0.020 


-0.030 ±0.022 


-0.041 ±0.022 


90 


0.643 ± 0.004 


0.647 ±0.004 


0.645 ± 0.004 


h' 


-0.055 ±0.010 


-0.039 ±0.012 


-0.047 ±0.012 



Table 1: Parameter best fit from two independent analyses (statistical error only) 



6 Other systematic uncertainties on the best fit parameters 

In addition to the systematic uncertainties associated with differences of the two analyses, 
the following potential sources of systematic errors have been considered (see Table: Ej) 

6. 1 Variation of the trigger efficiency over the Mqq fit interval 

The trigger efficiency has been measured using a sample of "minimum bias" events 
recorded continuously by a trigger requiring only the presence of a signal in at least 
two quadrants of the neutral hodoscope (during data taking the rate of this trigger was 
downscaled by a large factor). Within statistical errors the dependence of the trigger 
efficiency on Mqq is found to be consistent with the constant value etr = 0.928 ± 0.001 



for (2mo)^ < Mqq < 0.097 (GeV/c^)^. An equally good fit to the trigger efficiency is 
obtained using a 3^''^ degree polynomial. Varying the polynomial coefficients, so that the 
increases by an amount corresponding to ±lcr, (oq — 02)771+, a2m+, Qq and h' change as 
shown in Table El 

6.2 Dependence on the upper edge of the fit interval 

The upper edge of the Mqq fit interval has been varied from 0.094 to 0.107 (GeV/c^)^, 
resulting in variations of the best fit parameters with respect to the default upper bound 
= 0.097 (GeV/c2)2 as listed in tabled 

6. 3 Dependence on the position of the decay vertex 

As an additional check of the acceptance calculation, the K"^ 7r^7r°7r° events have 
been subdivided into two independent samples with the distance D of the reconstructed 
decay vertex from the LKr in the intervals 48 < D < 88 m, and 88 < -D < 136 
m, respectively. The best fit parameter values obtained from separate fits to the the two 
samples agree within statistics, providing no evidence for a possible systematic uncertainty 
associated with the position of the kaon decay vertex. 

6.4 Dependence on the charge sign 

The K"^ 7r±7r°7r° events consist of 1.470 x 10^ K+ and 0.817 x 10^ R-. Separate fits to 
these two samples give statistically consistent values for all best fit parameters. The two 
values of the slope parameter are go = 0.638 ± 0.005 for and go = 0.653 ± 0.006 for 
K~ , which disagree by ~ 1.9a. We take one half of their difference (0.008) as a systematic 
uncertainty on the value of go obtained by the fit to the full sample. 

6.5 Dependence on the distance between the vr^ track and the nearest photon 

The TT^ interaction in LKr may produce multiple energy clusters which are located, in 
general, near the impact point of the vr^ track and in some cases may be identified as 
photons. In order to study the effect of these fake photons on the best fit parameters 
we have repeated the analysis by varying the cut on the minimum distance d between 
each photon and the track impact point on LKr (both analyses A and B require d > 15 
cm). Varying d between 10 and 25 cm changes (ao — a2)m+ by ±0.002, while leaving 
the other parameters unchanged. We take this variation as a systematic uncertainty on 
(ao — a2)m+. 

6.6 Effect of LKr resolution and non-linear response at low photon energies 

The effect of possible uncertainties in the parameters describing the LKr energy resolu- 
tion has been simulated by adding a Gaussian noise with r.m.s. value of 0.06 GeV to 
the measured photon energies. An additional uncertainty may arise form the correction 
applied to the measured photon energies to account for the LKr non-linear response at low 
photon energies (typically < 2% at 3 GeV and becoming negligible above 10 GeV). The 
parameters describing this correction have been varied within limits chosen so that the 
measured 7r° mass for symmetric photon pairs does not depend on the 7r° energy. Varying 
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Figure 7: A = (data - fit) /data versus Mqq for various theoretical models: a) using the simple charge- 
exchange model of ref. J2J'; b) fit to the rescattering model of ref. '21; c) fit to the model of ref. 
including pionium formation; d) fit to the model of ref. 0/ excluding a 7 bin wide interval centred at 
Mqq — (2m+)^. The two vertical dotted lines in d) show the interval excluded form the fit. The point 

Mqq = (2™+)^ is indicated by the arrow. 



both the LKr resolution and non-hnearity correction parameters according to these pro- 
cedures changes (oq — 02)^2+ by ±0.001, while leaving the other parameters unchanged. 
We take this variation as an additional systematic uncertainty on (cq — 02)?^+ 

Table El lists all the systematics uncertainties discussed above. These are added in 
quadrature to obtain the total experimental systematic error on the values of the best fit 
parameters. 
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0.001 
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±0.009 
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0.009 


0.003 
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Table 2: Systematic uncertainties 



6.7 "External" uncertainties 

A crucial parameter in the model of refs. P [2] is the ratio R = A+^^/A+oo between 
the weak amplitudes of n'^n^'K^ and tt'^tt^tt^ decay. The value extracted 

from the measured decay branching ratio jj5j is i? = 1.972 ± 0.023. Varying R withing its 
error changes (oq — 02)?^+ by ±0.003, while leaving the other parameters unchanged. An 
additional theoretical error of ±5% on (ao — a2)m+, or ±0.013 is estimated in ref. |2j as 
the result of neglecting higher-order terms and radiative corrections in the rescattering 
model. These uncertainties have no significant effect on a2m_|_. 

Taking into account all systematic and external uncertainties we quote: 

(ao - a2)m+ = 0.268 ± O.OlO(stat) ± 0.004(syst) ± 0.013(ext) (3) 
a2m+ = -0.041 ± 0.022(stat) ± 0.014(syst) (4) 

The two statistical errors from the fit are strongly correlated, with a correlation coefficient 
of -0.858. We note that this analysis offers the first direct determination of a2, though 
not as precise as that of ao — a2. 

Preliminary results obtained under the assumption of exact isospin symmetry have 
been reported earlier |llj . 

7 Fit using the correlation between ao and a2 predicted by chiral symmetry 

It has been shown that analyticity and chiral symmetry provide a constraint between ao 
and a2 [T^ : 

a2m+ = (-0.0444±0.0008)±0.236(aom+ -0.22) -0.61(aom+ -0.22)2 -9.9(aom+ -0.22)^ 
Using this constraint in the fit to the rescattering model of ref. |2^ we obtain 

aom+ = 0.220 ± 0.006(stat) ± 0.004(syst) ± O.Oll(ext) (5) 
which corresponds to 

(ao - a2)m+ = 0.264 ± 0.006(stat) ± 0.004(syst) ± 0.013(ext.) (6) 



8 Summary and conclusions 



The TT^vr^ invariant mass (Mqo) distribution measured from a sample of 2.287 x 10^ — 
n'^n^'K^ fully reconstructed decays collected by the NA48/2 experiment at the CERN SPS 
shows an anomaly at Mqo = 2m^. This anomaly has been observed for the first time in 
this experiment thanks to the large statistical sample and the excellent Mqo resolution. 
It can be described by a rescattering model jT] |2I dominated by the contribution from 
the decay — *■ 7r^7r"'"7r~ through the charge-exchange reaction vr+vr" — > vr^vr^. These 
data have been used, therefore, to determine the difference — 02 between the J = 
and 1 = 2 S-wave tttt scattering lengths. Our result (see eq. Ej) is in very good agreement 
with theoretical calculations performed in the framework of Chiral Perturbation Theory 
(ChPT) |13j, which predict (ao — 02)^+ = 0.265±0.004. A different theoretical calculation 
based on a direct analysis of tttt scattering data without using chiral symmetry leads 
to a somewhat different value with a larger uncertainty, (oq — 02)^+ = 0.278 ± 0.016, 
which also agrees with our result. 

Previous determination of the vrvr scattering lengths have relied on a variety of meth- 
ods, such as the measurement of iT^7i~e~^h'e decay jlH], also being studied by the 
NA48/2 collaboration, or the measurement of the lifetime of the tt+tt" atom [l^j. Our 
value of Oo (see eq. El) is in good agreement with the result of experiment 865 at BNL |15j . 
aom+ = 0.216 ± 0.013(stat) ± 0.002(syst) ± 0.002(theor.), also obtained using constraints 
based on analyticity and chiral symmetry. Our value of oq — CI2 is also in good agreement 
with the first measurement of the lifetime of the vr+vr" atom ^H], which corresponds to 
|ao — a2\fn+ = 0.264l'lo;o2o (it should be noted that the latter result provides only a deter- 
mination of I ao — 02], while our measurement of —> n^n^n^ decays is also sensitive to 
its sign). 

To conclude, the study of a large sample of — > tt^tt^tt^ decays with excellent 
resolution on the Tr^Tr*^ invariant mass has provided a novel, precise determination of aQ—a2, 
independent of other methods and with different systematics uncertainties. In the near 
future the expected increase of the event sample by about a factor of 5 from the analysis of 
all the 2003-2004 data will further reduce the statistical error of our measurement. To be 
useful, this will require an improvement of the rescattering model to include higher-order 
terms and also radiative corrections. 
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